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We discuss the flavor dependence of the pressure and critical temperature cal- 
culated in QCD with 2, 2+1 and 3 flavors using improved gauge and staggered 
fermion actions on lattices with temporal extent N T = 4. For T J> 2T C we find 
that bulk thermodynamics of QCD with 2 light and a heavier strange quark is 
well described by 3-flavor QCD while the transition temperature is closer to that 
of 2-flavor QCD. Furthermore, we present evidence that the chiral critical point of 
3-flavor QCD, i.e. the second order endpoint of the line of first order chiral phase 
transitions, belongs to the universality class of the 3d Ising model. 



1 Introduction 

The existence of a finite temperature phase transition in strongly interacting 
matter is one of the most exciting non-perturbative features of QCD. De- 
termining the equation of state and the transition temperature is one of the 
basic goals in lattice studies of finite-T QCD. Studies of the transition which 
have been performed during the last years have shown that the details of the 
transition strongly depend on the number of quark flavors (n/) as well as the 
value of e.g. the pseudo-scalar meson mass (mps) which is controlled through 
variation of the bare quark masses (m q ). Furthermore, it became evident 
that the moderate values of the lattice spacing (a ~ 0.25fm) used in finite-T 
calculations with dynamical fermions lead to sizeable cut-off effects. Different 
discretization schemes used in the fermion sector, e.g. the standard stagger 
and Wilson fermion actions, lead to significantly different results for T c 
Calculations with improved actions thus seem to be mandatory to perform 
quantitative studies of the QCD equation of state at high temperature and to 
determine accurately the value of T c and its dependence on n/ and tops . We 
will report here on results for rif = 2 and 3 obtained with a Symanzik im- 
proved gauge action and an improved staggered fermion action, the p4-action 
with fat 1-link termscl. In addition we present results from a calculation with 
two light and a heavier strange quark. Details on the improved action as well 
as the algorithm used in our simulations are given in Ref. 2. 
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ERBFMRX-CT-970122 and the DFG grant Ka 1198/4-1. 
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Figure 1. The pressure in units of T 4 (left) and psB (right) for 2, 2+1 and 3-flavor QCD 
obtained from calculations with the p4 action on lattices with temporal extent N T = 4. 
The light quark mass used in all cases is m q /T = 0.4 and the heavier quark mass used in 
the (2+l)-flavor case is m s /T = 1. 



2 Flavor dependence of the QCD equation of state 

In the high (infinite) temperature limit the QCD equation of state is expected 
to approach that of an ideal quark-gluon gas, i.e. bulk thermodynamic observ- 
ables like energy density and pressure will reflect the number of light degrees 
of freedom, 

esB /r 4 = 3psb/t 4 = (s + fs/)^ , (i) 

where g f = Ei=«,rf,.. 9(m/T) with 

g(m/T) = —r dxxJx 2 - (m/T) 2 ln(l + e~ x ) , (2) 

?7T 4 Jm/T 

counts the effective number of degrees of freedom of a massive Fermi gas. 
For a gas of massless quarks one has, of course, gf = rif. The effective 
number of degrees of freedom, g(m/T), rapidly approaches unity for quark 
masses smaller than T. For instance, one has g(l) — 0.8275 and g(0.4) = 
0.9672. These jjiumbers correspond to the bare quark mass values used in 
our simulationso on lattices of size 16 3 x 4. Results obtained for the pressure 
are shown in Fig. [I]. In the high-T phase p/T A clearly shows the expected 
flavor dependence. We note that lattice calculations are performed at fixed 
rriq/T . Using instead a fixed physical quark mass value, e.g. m s ~ T c , would 
require to reduce m q /T in a simulation as the temperature is increased. From 
Fig. |l| we thus conclude that already at T ~ 2T C the pressure calculated in 
(2+l)-flavor QCD with a fixed strange quark mass will be similar to that of 
massless 3-flavor QCD. 
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Figure 2. The critical temperature in units of y/a (left) and the vector meson mass (right) 
versus (mpg/my) 2 . Shown are results for 2, (2+1) and 3-flavor QCD obtained from cal- 
culations with the p4 action on lattices with temporal extent N T = 4. For rif = 2 we also 
show results obtained by using unimproved gauge and staggered fermion actions. 



3 Flavor and quark mass dependence of T c 

The transition temperature in QCD with dynamical fermions has been found 
to be significantly smaller than in the pure gauge sector. This is in accordance 
with intuitive pictures of the phase transition based e.g. on the thermody- 
namics of bag or percolation models. With decreasing m q the hadrons become 
lighter and it becomes easier to build up a sufficiently high particle density 
that can trigger a phase transition. For the same reason such models also 
suggest that T c becomes smaller when n/ and in turn the number of light 
pseudo-scalar mesons increases. This qualitative picture is confirmed by the 
lattice results presented below. 

We have determined the pseudo-critical couplings, j3 c {m q ), for the tran- 
sition to the high temperature phase on lattices of size 16 3 x 4. For 2 and 
3-flavor QCD calculations have been performed in a wide range of quark 
masses, 0.025 < m q < 1.0. The smallest quark mass corresponds to a pseudo- 
scalar meson mass raps — 350 MeV. In order to set a scale for the transition 
temperature we calculated the light meson spectrum and the string tension^] 
at f3 c (m q ) on lattices of size 16 4 . The resulting pseudo-critical temperatures 
are shown in Fig. ^|. 

We note that T c /y/a and T c /mv do show a consistent flavor dependence 
of T c . In a wide range of quark mass values, T c in 3-flavor QCD is about 

a We note that the heavy quark potential is no longer strictly confining in the presence 
of dynamical quarks. The definition of the string tension is based on potentials extracted 
from Wilson loops which have been found to lead to a linear rising potential at least up to 
distances r ~ 2fm. 
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10% smaller than in 2-flavor QCD. The dependence on m q is, however, quite 
different when using the vector meson mass rather than the string tension to 
set the scale for T c . While T c /yfa does show the expected rise with increasing 
values of m q the contrary is the case for T c jmy. Of course, this does not 
come as a surprise. The vector meson mass used in Fig. || to set the scale 
is itself strongly quark mass dependent, my = m p + c m q . Its mass thus 
is significantly larger than the vector meson mass, m p , in the chiral limit. 
We stress this well known fact here because it makes evident that one has to 
be careful when discussing the dependence of T c on parameters of the QCD 
Lagrangian, e.g. nj and m q . One has to make sure that the observable used 
to set the scale for T c itself is not or at most only weakly dependent on the 
external parameters. The fact, that the hadron spectrum as well as the string 
tension calculated in quenched (m q — > oo) QCD0 are in reasonable agreement 
with experiment and phenomenology and show systematic deviations only on 
the 10% level indicates that the corresponding partially quenched observables 
are suitable also for setting the scale in the presence of dynamical quarks. 

4 T c in the chiral limit 

The chiral phase transition in 3-flavor QCD is known to be first order whereas 
it is most likely a continuous transition for nt = 2. This also implies that the 
dependence of the (pseudo)-critical temperature on the quark mass will differ 
in both cases. Asymptotically, i.e. to leading order in the quark mass one 
expects to find, 



with 1/(35 — 0.55 if the 2-flavor transition indeed belongs to the universality 
class of 3d, 0(4) symmetric spin models. Our estimates of T c in the chiral 
limit are based on the data shown in Fig. 2. In the case of «/ = 3 we have 
extrapolated T c /y/a and T c /mv using an ansatz quadratic in mps/my. In 
addition we have extrapolated my calculated at (3 c (m q ) for m q = 0.025 and 
0.05 linearly in m q to the critical point in the chiral limit, f3 c (0) — 3.258(4). 
The extrapolation to the chiral limit is less straightforward in the case of 
rif = 2. The data shown in Fig. 2 indicate that the quark mass dependence 
for rif — 2 and 3 is quite similar. This suggests that sub- leading corrections, 

E. g. the ratio y/a/m p calculated in the limit of vanishing valence quark mass does show 
little dependence on the dynamical quark mass used to generate gauge field configurations. 




, n f = 2 
, n f > 3 



(3) 
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quadratic in mps I my , should be taken into account in addition to the leading 
behavior expected from universality. In our extrapolations for nj = 2 we thus 
also add a term quadratic in mps/my to the leading term given in Eq. 3. 
From these fits we estimate 

T c _ J 0.225 ±0.010 ,n f = 2 

0.20 ±0.01 , n/ = 3 ' (4) 

which corresponds to T c = (173±8) MeV and (154±8) MeV for n f = 2 and 3, 
respectively. For n/ = 2 this agnees well with results obtained in a calculation 
with improved Wilson fermionsu. We stress, however, that the errors given 
here as well as in Ref. 3 are statistical only. Systematic errors resulting from 
remaining cut-off effects and from the ansatz used for extrapolating to the 
chiral limit are expected to be of similar magnitude. In order to control these 
errors calculations on lattices with larger temporal extent are still needed. 

Physically most relevant is a determination of the transition temperature 
of QCD with 2 light quarks (m u ^ ~ 0) and a heavier (strange) quark with 
m s ~ T c . The result from our calculation with m u .d/T = 0.4 and m s /T = 1 
is also shown in Fig. 2. Although this analysis of (2+l)-flavor QCD has 
not yet been performed with a sufficiently light light quark sector the result 
obtained for the pseudo-critical temperature does suggest that also in the case 
of realistic light quark masses T c will be close to that of 2-flavor QCD. 

5 Universality at the chiral critical point of 3-flavor QCD 

In the chiral limit of 3-flavor QCD the phase transition is first orderO. It thus 
will persist to be first order for m q > up to a critical value of the quark 
mass, m c . At this chiral critical point the transition will be second order. It 
has been conjectured that this critical point belongs to the universality class 
of the 3d Ising modelQ. From a simulation with standard, i.e. unimproved, 
gauge and staggered fermion actions we find support for this conjecture. 

As discussed in the context of the electroweak transition the analysis of 
the critical behavior at the 2 nd order endpoint of a line of 1 st order phase 
transitions requires the correct identification of energy-like and ordering-field 
(magnetic) directional A general approach to determine the corresponding 
operators has been discussed in the context of the 3d, 3-state Potts model! 
The ordering field operator at the chiral critical point can be constructed from 
a linear combination of the gluonic action Sq and the chiral condensate i^tj), 

E = S G + r^tp , M = tpip + sS G . (5) 

Here the mixing parameter r is determined from the m^-dependence of the line 
of first order transitions, r = (d/3/dmg)^ point and s by demanding (E ■ M) — 
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Figure 3. Contour plot for the joint probability distribution of Sc and 1(11(1 (left) as well as 
E and M (middle) for 3-fiavor QCD and the 3d, 3-state Potts model (right). The QCD 
contour plots are based on calculations performed on a 16 3 X 4 lattice at = 5.1499 with 
m q = 0.035. The mixing parameters have been fixed to r = 0.55 and s = 0. Parameters 
for the Potts model are given in Ref. 7. 



0. In fact, unlike for energy-like observables, e.g. critical amplitudes that 
involve the thermal exponent yt , the universal properties of observables related 
to M do not depend on the correct choice of the mixing parameter s as 
long as the magnetic exponent yn is larger than yt. The joint probability 
distribution for E and M characterizes the symmetry at the critical endpoint 
and its universality class. Contour plots for the E-M distributions at the 
critical endpoints of 3-flavor QCD and the 3d, 3-state Potts model as well 
as the corresponding plot for the Sc-i^tp) distribution are shown in Fig. 3. 
This shows that also in the QCD case a proper definition of energy-like and 
ordering-field operators is needed to reveal the symmetry properties of the 
chiral critical point. The joint distributions of E and M suggest the universal 
structure of the E-M probability distribution of the 3-d Ising modelS, although 
it is apparent that an analysis of 3-flavor QCD on larger lattices is needed to 
clearly see the "two wings" characteristic for the 3d Ising distribution. 
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